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Abstract. Unitary matrices which are zero below the secondary diagonal 
(Hessenberg unitary matrices) have many uses in analysis. Given a set of 
needed conditions on a unitary matrix, this algorithm will give the sparsest 
unitary matrix. We give an algorithm for constructing all real Hessenberg 
unitary matrices. The n X n unitary matrices given by the algorithm have 
n — 1 variables which can be chosen to give additional properties needed for a 
particular application. 



1. Introduction 

Hessenberg matrices have may uses in both pure an applied mathematics. It 
is a well used fact that every matrix is unitarily equivalent to such a matrix and 
there are polynomial time algorithms for doing this. In numerical linear algebra, 
they are used to speed-up eigenvalue computations. Also, real Hessenberg unitary 
matrices have various uses in analysis. They also have the advantage of being quite 
sparse. Here we give an algorithm which constructs all such real unitary matrices. 
The algorithm starts with 2x2 matrices and for each new dimension, we drop the 
top row of the previous dimension, introduce a new variable and then construct 
two new rows with the new variable added and finally use the rest of the rows of 
the previous case by adding the first column of zeroes. The n x n unitary matrices 
have n — 1 variables which can be chosen to give additional properties needed for a 
particular application. 

This algorithm was developed at the request of the people working in the Frame 
Research Center (www.framerc.org) for their work in frame theory where they 
needed the sparsest unitary matrices. 

Remark 1.1. Since we can multiply rows and columns of a unitary by —1 and 
still have a unitary, and if we permute rows or columns while maintaining the 
Hessenberg property we still have a unitary, we will construct all real triangular 
unitary matrices up to such changes in sign and permutations of rows or columns. 



2. The Construction 

We will give an algorithm which constructs all real Hessenberg unitary matrices. 
So we can better see how the algorithm works, we will start by giving some small 
examples. 
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2.1. 2x2 Matrices. This is a unique class and is given by: 

Here we must have < z\ < 1. 

2.2. 3x3 Matrices. This class is given by: 

V 1 - z 2 y/(l ~ gljfg V-gl^2 

~V^2 - Z 2 )(l - Zi) - Z 2 ) zi 

-yfZl V 1 - z l 

Here, < zi, z 2 < 1. 

Proposition 2.1. T/iis c/ass of matrices consists of unitary matrices. 
Proof. The proof is done by cases: 

Case 1: The rows have norm 1. 

The square sum of the elements of row 1 is 

(1 - z 2 ) + (1 - zi)z 2 + ZlZ 2 = 1. 
The square sum of the elements of row 2 is 

z-i + (1 - z 2 )(l - zx) + (1 - z 2 )z 1 = z 2 + (l- z 2 ) = 1. 
It is clear that row 3 is norm 1. 

Because of the symmetry of the matrix (I.e. switching Z\ and z 2 in the calcula- 
tion), it is clear that the row sums and column sums are equal. 

Case 2: The rows are orthogonal. 

The inner product of rows 1 and 2 is 

-\J{\ - Z 2 )z 2 + (1 - Zl)y/(1 - Z 2 )Z 2 + 21^/(1 - Z 2 )Z 2 = 

y/(l-z 2 )z 2 [-\ + (1 - Zl ) + zi] - 0. 
The inner product of rows 1 and 3 is 

-\J (1 - zi)ziz 2 + - zi)ziz 2 = 0. 
The inner product of rows 2 and 3 is 

-y/(l - z 2 )(l - zi)zi + y/(l - z 2 )(l - Zl ) Zl = 0. 

By symmetry again, the columns are orthogonal. □ 

We give some examples of how to use this: 
Case 1: If z\ = z 2 = we get 

"1 0" 

1 
1 
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Case 2: If we let z\ = z 2 = 1 we get 

"0 1 

-10 
0-10 

Case 3: If wc let z\ = 1 and Z2 — we get 

"l 0" 

1 
0-10 

Case 4: If we let z\ = \ and z 2 = 1 we get 
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Case 5: If z\ = 




Theorem 2.2. This algorithm gives all 3 x 3 real Hessenberg unitary matrices. 

Proof. We will just outline this. Given a unitary 

an 0,12 ai3 
a 2 i a 22 «23 
_ a 32 a 33 _ 

If one of a 32 or a 33 is zero, it is easily checked that the matrix is a permutation of 
the rows or columns of the identity multiplied by —1 if necessary and this matrix 
easily arises from our algorithm. Otherwise, we may assume a 32 is negative and 
let our z\ = a\ 2 to get the last row correct given that it square sums to 1. Now, 
the last row must be orthogonal to the rows above it so (a 2 2 ? <223) = c(ai 2 ,ai 3 ) are 
the unique vectors (up to length) in K 2 which are orthogonal to row 3. Since z\ is 
given, we can find the unique z 2 which makes the last two terms of our first and 
second rows equal to the given one. The first column is now uniquely determined 
and must be of our form. □ 



2.3. 4x4 Matrices. We construct a unitary matrix by: 



VI - z 3 y/(l - Z 2 )Z 3 - Zl)z 2 Z 3 y/ZlZ 2 Z 3 

~\/Z3 Vi 1 - z 3)(l - Z 2 ) y/( l - Z 3 )(l - Zi)z 2 - Z 3 )z~ 2 

-y/zi ^(1-Z 2 )(l- Zl ) y/( l-Z 2 ) Zl 

-JI{ VT^~zi 



Wc leave it to the reader to check that this constructs all 4 x 4 real Hessenberg 
unitary matrices. 
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2.4. 5x5 Matrices: The Algorithm. Given our general n x n real Hcsscnbcrg 
unitary matrix with n — 1 variables in it, we introduce a new multiplier y/z n+ i 
on the second row of the n x n case and the first entry in the second row of the 
(n + 1) x (n + 1) matrix is — ^/2„+i. We also need the sum of the squares of the 
entries of the two new rows we are adding at the top to equal the squares of the 
entries of the top row of the n x n matrix. After that, we add the remaining rows 
of the previous case with a column of zeroes in front. 

So, to pass from the 4x4 case to the 5x5 case we create two new rows: 



V 1 - z i y/Q- - Z 3 )z 4 - Z 2 )z 3 Z 4 - 2i)2 2 2 3 2 4 y / Z 1 Z 2 Z 3 Z 4 

~\fzi Vi 1 - Z 4)(l - Z 3 ) - Z 4 )(l - Z 2 )z 3 - 2 4 )(1 - Zi)z 2 Z 3 Z A )ziZ 2 Z 3 

Combining this with the previous case we get: 



VI - z 4 - 23)24 V(l - z 2 )z 3 z 4 V(l - z 1 )z 2 z 3 z 4 ^z x z 2 z 3 z 4 ' 

V(l - Z 4 )(l - Z 3 ) Z 4 )(l ~ Z 2 )z 3 Z 4 )(l - Zi)z 2 Z 3 y/ (1 - Z 4 )Z\Z 2 Z 3 



Z 4 _ _ _ 

-ypti ^{l- Z 3 ){1- Z 2 ) y/{ \- Z 3 ){l-^y 2 - Z 3 ) Zl Z 2 

-y/zi ^(l-z 2 )(l- Zl ) ^( l-z 2 ) Zl 

-V^T V^zi 

Theorem 2.3. The matrix given above is a unitary matrix. 

Proof. We will outline the proof. Since the last three rows come from the previous 
unitary matrix, we just need to check that the first two rows are orthogonal to all 
the others, that they are norm 1 and the column vectors have norm 1. 

The two new rows are norm 1: 

The norm of row 1 is: 

1 - z 4 + (1 - z 3 )z 4 + (1 - z 2 )z 3 z 4 + (1 - 21)2:22:32:4 + 212:22:32:4 = 
1 - 2324 + (1 - 22)2324 + Z 2 Z 3 Z 4 = 1 - 2 3 2 4 + 2 3 2 4 = 1. 

The norm of row 2 is: 

2 4 + (1 - 2 4 )(1 - 2 3 ) + (1 - 2 4 )(1 - 22)23 + (1 - 2 4 )(1 - 21)2223 + (1 - 2 4 )2l2 2 2 3 = 
2 4 + (1 - 2 4 ) [(1 - 2 3 ) + (1 - 22)23 + (1 - 2!)2 2 2 3 + 2l2 2 2 3 ] = 
2 4 + (1 - 2 4 ) [1 - 2 2 2 3 + 2 2 2 3 ] = 1. 

The column vectors are norm 1: 

Column 1 square sums to (1 — 24) + 24 = 1. Column 2 square sums to 

(1 - 23)24 + (1 - 2 4 )(1 - 2 3 ) + 2 3 = 1 - 2 3 + 2 3 = 1. 

Column 3 square sums to 

(1 - 22)2324 + (1 - 2 4 )(1 - 22)23 + (1 - 2 3 )(1 - 2 2 ) + 2 2 = 
(1 - 2 2 ) [2324 + (1 - 24)23 + (1 - 2 3 )] + 2 2 = (1 - 2 2 ) [1] + 2 2 = 1. 

Column 4 square sums to 

(1 - 2l)2 2 2 3 2 4 + (1 - 2 4 )(1 - 21)2223 + (1 - 2 3 )(1 - 2l)22 + (1 - 2 2 )(1 - Z\) + Z\ = 
(1 - 2l) [2 2 2 3 2 4 + (1 - 24)2223 + (1 - 23)22 + (1 - Z 2 )\ + 2l = 
(1 - 2l) [2 2 2 3 + 1 - 2 2 2 3 ] + 2l = 1. 
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Column 5 square sums to 

z x z 2 z 3 z 4 + (1 - z 4 )ziz 2 z 3 + (1 - z 3 )ziz 2 + (1 - z 2 )zi + (1 - Zl) = 

Z1Z2Z3 + ZxZ 2 - ZlZ 2 Z 3 + X\ - ZlZ 2 + 1 - Zl = 1. 

The rows are orthogonal: The inner product of rows 1 and 2 is: 

-V/(Z 4 (1 - 24) + y/(l - Z 3 ) 2 Z 4 (1 - Zi ) + yJ(l- Z 2 ) 2 Z 2 Z 4 (1 - 24) + 

y^i - zi) 2 zf zf z 4 (i - z 4 ) + y / z 2 z|z|z 4 (i - Z 4 ) = 

-v/z 4 (l - z 4 ) + (1 - z 3 ) v /z 4 (l - z 4 ) + (1 - z 2 )z 3V /z 4 (l - z 4 )+ 
(1 - zi)z 2 z 3 ^z 4 (l - z 4 ) + ziz 2 z 3 - v /z 4 (l - z 4 ) = 
y/z±(l - z 4 )[-l + (1 - z 3 ) + (1 - z 2 )z 3 + (1 - zi)z 2 z 3 + Z1Z2Z3] = 
\j z 4 (l - z 4 )[-l + (1 - z 3 ) + z 3 - z 2 z 3 + z 2 z 3 - ZlZ 2 Z 3 + ZlZ 2 Z 3 ] = 0. 

□ 

The inner product of rows 1 and 3 and taking our squares as above at the same 
time is 

-^3(1 - z 3 )z 4 +(l-z 2 )v / z 3 (l - z 3 )z 4 +(1-zi)z 2 ^:e 3 (1 - z 3 )z 4 +ziz 2 ^z 3 (l - Z 3 )Zi = 

^z 3 (l-z 3 )z 4 [-l + (1 - z 2 ) + (1 - zi)z 2 + ziz 2 ] = 
The inner product of rows 1 and 4 is 

-\Jz 2 (l - z 2 )z 3 z 4 + (1 - Z X )\J (1 - z 2 )z 2 z 3 z 4 + Z 1 y / (1 - z 2 )z 2 z 3 z 4 = 

v /z 2 (l-z 2 )z 3 z 4 [-l + (1 - Zl) + Zl] = 0. 
The inner product of rows 1 and 5 is 

-\/zi(l - zi)z 2 z 3 z 4 + \J Zl(l - zi)z 2 z 3 z 4 ] = 0. 
The inner product of rows 2 and 3 is 

-y/(l - z 4 )(l - z 3 )z 3 + (1 - z 2 )v/(l - z 4 )(l - z 3 )z 3 + 

(1 - z 1 )z 2 y / (l - z 4 )(l - z 3 )z 3 + ziz 2V / (1 - z 4 )(l - z 3 )z 3 = 

v/(l-z 4 )(l-z 3 )z 3 )[-l + (1 - z 2 ) + (1 - zi)z 2 + ZlZ 2 ] = 0. 
The inner product of rows 2 and 4 is 

- z 4 )z 2 (l - z 2 )z 3 +(l-zi)\/(l - z 4 )(l - z 2 )z 2 z 3 +ziv/(l - z 4 )(l - z 2 )z 2 z 3 = 0. 
The inner product of rows 2 and 5 is 

- z 4 )zi(l - zi)z 2 z 3 + \/(l - z 4 )zi(l - zi)z 2 z 3 = 0. 
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3. The General Case 

For the general case, we take the n x n case and delete the first row and add 
a first column of zeroes. Then we add two new rows at the top of this matrix by 
choosing a new variable z n and making row 1 as: 

Vl - z n - z n ^ijz n z n -2)z n -iz n - z n - 3 )z n -iz n - 2 ■■■ 

— Z 1 )z 2 Z 3 • • • z n ■ ■ ■ \fZ\Z 2 ■ ■ ■ Z n 

and making row 2 as 

-\f&k V(l - Z n )(l - Z n -{) - z n )(l - z n - 2 )z„-i - z n )(l - z n ^ 3 )z n -iz n - 2 ■ ■ ■ 

- z n )ziz 2 z 3 . . . Z n -i. 

Remark 3.1. It takes a significant amount of effort to show that we can produce 
all the required Hessenberg unitary matrices and it does not seem to be important 
enough to justify the effort since our intention is not to publish this but just post 
it on the arXiv so it is available to researchers. So we will not address this here. 
Basically, it can be done by induction on n and a case analysis of the placement of 
zeroes. 
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